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Abstract 

The energy dependence of charged- hadron production in relativistic heavy- 
ion collisions is investigated in a nonequilibrium-statistical relativistic diffu- 
sion model (RDM) with three sources. Theoretical pseudorapidity distribu- 
tions are compared with Au + Au data at RHIC energies of v /s/vjv = 0.13 
and 0.2 TeV, and computed for Pb + Pb central collisions at LHC energies 
of 2.76 and 5.52 TeV. The nearly equilibrated source at midrapidity arising 
from gluon-gluon collisions becomes the major origin of particle production 
at LHC energies. The midrapidity dip is determined by the interplay of the 
three sources. 
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1. Introduction 

With the advent of first results from heavy-ion collisions at LHC energies 
of a/s/vtv = 2.76 TeV in central Pb + Pb collisions [U 12], a new perspective 
on this area of research opens up. The strong gluon field that is present at 
these high energies determines the dynamics of the collision and the details of 
particle production even more decisively than in Au + Au collisions at RHIC 
energies of 0.13 and 0.2 TeV, where quark-gluon interactions are still more 
important in the particle production process than gluon-gluon collisions. The 
first and simplest observable to be determined experimentally is the charged- 
particle multiplicity density at mid-rapidity in central Pb + Pb. There are 
many theoretical models predicting this value with varying accuracy (see (3j 
H] at the maximum LHC energy of 5.52 TeV, and [1] at 2.76 TeV). However, 
the experimental ALICE result of 1601 ±60 [2] at 2.76 TeV is obtained from 
a straightforward extrapolation of the midrapidity values at RHIC energies 
with log(y'sjvjv)- 

More specific information can be expected from the detailed shape of 
the pseudorapidity distribution of produced charged hadrons at 77— values 
further away from midrapidity, which will be available experimentally in 
the near future. The decomposition of the distribution function (dN/drj)(rj) 
from the underlying physical ingredients such as quark-gluon vs. gluon-gluon 
interactions will be of particular interest. 

In this Letter an analytically soluble nonequilibrium-statistical RDM- 
model [5j |6] that successfully describes pseudorapidity distributions for pro- 
duced hadrons at RHIC energies is used to predict these distribution func- 
tions at LHC energies. The model relies on three sources for charged-hadron 
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production, with the midrapidity source associated with gluon-gluon col- 
lisions, and two forward-centered fragmentation sources arising essentially 
from valence quark - gluon interactions. 

It has been shown in [3 El E] within the relativistic diffusion model (RDM) 
that at RHIC energies of 0.13 TeV (0.2 TeV) the midrapidity source gener- 
ates about 13 % (26 %) of the produced particles in a 0-6% central Au + Au 
collision, whereas the bulk of the particles is still produced in the two frag- 
mentation sources. At SPS, and low RHIC energies of 19.6 GeV the effect of 
the midrapidity source is negligible [9]. 

In the asymmetric d + Au system at 0.2 TeV there is also a sizeable 
midrapidity source containing 19 % of the produced particles for 0-20% cen- 
tral collisions [in]. Particle creation from a gluon-dominated midrapidity 
source, incoherently added to the sources related to the valence part of the 
nucleons, had also been proposed by Bialas and Czyz [TTJ. There exist also 
many other models which assume a central source such as the dual parton 
model [T2| 113] . or the quark-gluon string model [H]. The RDM provides an 
analytical framework to investigate the interplay of central and fragmenta- 
tion sources. 

For asymmetric systems, the central source is shifting in rapidity space 
with increasing centrality, whereas for symmetric systems it remains at midra- 
pidity < r] >= 0. The shape of the d iV7 (^-distributions at different cen- 
tralizes is very sensitive to the detailed balance of the underlying distribu- 
tion functions, and the excellent agreement with the d + Au PHOBOS-data 
[T5| [T6| ITT] at 0.2 TeV lends credibility to the three-sources model also for 
symmetric systems where the details of the distribution functions are less 
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specific. 

Within the RDM, I investigate in this Letter the energy dependence of 
the three sources for particle production in central collisions of symmetric 
systems, and provide predictions at LHC energies. The energy range consid- 
ered here for the three-sources model covers RHIC energies of ^snn = 0.13 
and 0.2 TeV in Au + Au collisions, the presently accessible LHC energy of 
2.76 TeV in Pb + Pb collisions, and the maximum LHC energy of 5.52 TeV. 

The model is considered in Sec. 2, the calculation of pseudorapidity 
distributions of charged hadrons at RHIC and LHC energies in Sec. 3, and 
conclusions are drawn in Sec. 4. 

2. Relativistic Diffusion Model 

In the Relativistic Diffusion Model, the rapidity distribution of produced 
particles emerges from an incoherent superposition of the beam-like fragmen- 
tation components at larger rapidities arising mostly from valence quark- 
gluon interactions, and a component centered at midrapidity that is essen- 
tially due to gluon-gluon collisions. All three distributions are broadened in 
rapidity space as a consequence of diffusion-like processes. 

The time evolution of the distribution functions is governed by a Fokker- 
Planck equation (FPE) in rapidity space [5j (and references therein) 



<9 r „ d 
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dt [R(y,t)r = ~-[J(y)[R(y,t)r\+—[D y .R(y,tW (1) 

with the rapidity y = 0.5 • ln((E + p)/(E — p)). The beam rapidity can 
also be written as ybeam = TVmax = ^^{y/sNN / rn p ) . The rapidity diffusion 
coefficient D y that contains the microscopic physics accounts for the broad- 
ening of the rapidity distributions. The drift J(y) determines the shift of the 
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mean rapidities towards the central value, and linear and nonlinear forms 
have been discussed [TH1 EE E] . 

The standard linear FPE corresponds to fi — v — 1 and a linear drift 
function 

J{y) = (y eq -y)/T y (2) 

with the rapidity relaxation time T y , and the equilibrium value y eq of the 
rapidity. This is the so-called Uhlenbeck-Ornstein [20] process, applied to 
the relativistic invariant rapidity for the three components Rk{y, t) (k— 1,2,3) 
of the distribution function in rapidity space 



d „ , , Id 
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dt R k (y,t) = --—^ Veq - y ) ■ R h {y,t)\ + — [D k y ■ Rk(y,t)\ . (3) 

Since the equation is linear, a superposition of the distribution functions 
[SJ EI] using the initial conditions Ri^iy, t = 0) = 8(y ± y ma x) with the 
absolute value of the beam rapidities y ma x, and R^{y,t = 0) = S(y — y eq ) 
yields the exact solution. In the solution, the mean values are obtained 
analytically from the moments equations as 

< Vifi(t) >= y eq [l - exp(-t/T y )) =f y max exp {-t/r y ) (4) 

for the sources (1) and (2) with the absolute value of the beam rapidity 
y ma x, and y eq for the local equilibrium source which is equal to zero only for 
symmetric systems. Hence, both mean values < y± t 2 > would attain y eq for 
t— > oo, whereas for short times they remain between beam and equilibrium 
values. The variances are 

*W*) = Dy^yl 1 ~ exp(-2t/r y )], (5) 



6 



and the corresponding FWHM- values are obtained from T = a/8 In 2 • a since 
the partial distribution functions are Gaussians in rapidity space (but not in 
pseudorapidity space). 

The midrapidity source has mean value zero and hence, comes close to 
thermal equilibrium with respect to the variable rapidity during the interac- 
tion time Ti n t. Note that the width approaches equilibrium twice as fast as the 
mean value. I use the notion R eq (y,t) for the associated partial distribution 
function in y-space, with charged particles, cf. Table [TJ Full equilibrium 
as determined by the temperature would be reached for T int /T y ^> 1. The 
fragmentation sources do not reach < y 12 >= during the interaction time 
and hence, remain far from thermal distributions in rapidity space, and do 
not fully equilibrate with the central source. 

3. Pseudorapidity distributions 

If particle identification is not available, one has to convert the results to 
pseudorapidity, rj = — ln[tan(#/2)] with the scattering angle 9. The conver- 
sion from y— to rj— space of the rapidity density 



The average mass < m > of produced charged hadrons in the central 
region is approximated by the pion mass since pions represent by far the 
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is performed through the Jacobian 



J(r), (m)/(pr» = cosh(r/)- 



[l + ((m)/(p T >) 2 + sinh 2 (77)] 



1/2 



(7) 



Table 1: Three-sources RDM-parameters for 0-6% Au + Aii at RHIC energies (upper two 
lines) and for 0-5% Pb + Pb at LHC energies (lower two lines). See Fig. [2] and text for the 
extrapolation of the time parameter Ti n t/T v to LHC energies. Widths and particle numbers 
denoted by * are extrapolated linearly with log(y / s]v]v). At RHIC energies the nonequilib- 
rium sources from quark-gluon interactions with particle content A^ 2 dominate. At LHC 
energies the local equilibrium source from gluon-gluon collisions with particle content N^i 
is the major origin of particle production at midrapidity. Experimental midrapidity values 
(last column) are from PHOBOS [H [23] for \rj\ < 1 at RHIC energies and from ALICE 
for |ry| < 0.5 at 2.76 TeV. 
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largest fraction of produced charged hadrons, in particular in the midrapidity 
source where the transformation has the biggest effect. 

The dependence on the mean transverse momentum < > is illustrated 
in Fig. [TJ Due to the Jacobian, the partial distribution functions differ from 
Gaussians. In the actual calculations, I use < >= 0.3 and 0.4 GeV at the 
respective RHIC energies of 0.13 and 0.2 TeV, and < pt >= 0.6 and 0.7 GeV 
at LHC energies of 2.76 and 5.52 TeV. The values at LHC energies should 
be updated once measured py-distributions become available. 
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Figure 2: Dependence of the diffusion-model parameters for heavy systems (central Au 
+ Au at RHIC energies, central Pb + Pb at LHC energies) on the center-of-mass energy 
\J s nn according to [S]: Quotient of interaction time and relaxation time for sinh- and 
exponential (dashed) extrapolation (upper frame, with rescaled absolute values); width 
of the peripheral sources including collective expansion (middle frame); effective width of 
the midrapidity source (lower frame). The results are for charged-hadron pseudorapidity 
distributions, with extrapolations to LHC energies. The dots refer to the fit values at 
RHIC energies of 19.6, 130 and 200 GeV. The time parameters used in the present work 
at LHC energies of 2.76 and 5.52 TeV have been averaged between the two analytical 
extrapolations. 
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Figure 3: Calculated pseudorapidity distributions of produced charged particles from Au 
+ Au collisions (bottom) at ^snn — 0.13 and 0.2 TeV for 0-6% central collisions in 
comparison with PHOBOS data [HI [53]. The analytical three-sources RDM-solutions are 
optimized in a fit to the data. Distribution functions for 0-5% central Pb + Pb collisions 
at LHC energies of 2.76 and 5.52 TeV are shown in the upper part of the figure, with 
the lower-energy result adjusted to the recent midrapidity ALICE data point [2]. Dotted 
curves are without the Jacobian transformation. The corresponding parameter values are 
given in Table [T] 
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Figure 4: Pseudorapidity distributions of charged hadrons in 0-5% central Pb + Pb colli- 
sions at LHC energies of ^snn = 2.76 and 5.52 TeV. The underlying theoretical distribu- 
tions are shown for 2.76 TeV. Their shapes are not significantly modified by the Jacobian. 
The size of the midrapidity dip is determined by the interplay of central (gluon-gluon, 
dashed; without Jacobian, dotted) and peripheral (valence quarks - gluon, dash-dotted) 
distribution functions. The midrapidity value is almost completely determined by particle 
production from gluon-gluon collisions at LHC energies. 
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The dependencies of the diffusion-model parameters on incident energy, 
mass and centrality at RHIC energies have been investigated for various 
systems in [HI El 121] ■ In particular, the centrality dependence seen in the 
RHIC data is exactly reproduced [El E]- The parameters are shown in Fig. [2] 
and Table [T] as functions of the cm. energy in central collisions of Au + Au, 
and in an extrapolation to Pb + Pb at LHC energies. The difference between 
these two systems is very small since the diffusion- mo del parameters scale 
with the extension of the system like A 1 / 3 , which differs only by a factor of 
1.02. 

The time parameter T int /r y is displayed as function of center-of-mass en- 
ergy in the upper frame of Fig. [2j with a functional dependence on the beam 
rapidity y beam and hence, on energy given by T int /r y oc y b eamN part / smh(y beam ) 
as motivated in [21], whereas the dashed curve assumes an exponential de- 
pendence that yields a broader distribution function, see Fig. 8 in [21] for a 
detailed comparison of the two limiting cases. At LHC energies of 2.76 and 
5.52 TeV for Pb + Pb I use in this prediction intermediate values between 
the two analytical extrapolations, see Table [T] 

The partial widths (FWHM) as functions of energy within the RHIC 
range for Au + Au are displayed in the middle and lower frames of Fig. [2] 
for both fragmentation and midrapidity sources. Here the widths are ef- 
fective values: beyond the statistical widths that can be calculated from a 
dissipation-fluctuation theorem [22] within the RDM, they include the effect 
of collective expansion. The values at RHIC energies are resulting from a 
X 2 -minimization with respect to the data that corresponds to the time evo- 
lution up to Tint'. The integration is stopped at the optimum values of T int /T y , 
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Ti^eq, and N^; the explicit value of r int is not needed. 

The normalization is given by the total number of produced charged 
hadrons that is taken from experiment if available, or extrapolated in case of 
predictions at higher energies. Hence, the model contains five parameters for 
symmetric systems, and six parameters for asymmetric systems. It provides 
an analytical framework to calculate the distribution function, and to draw 
physical conclusions. 

The charged-particle distribution in rapidity space is obtained as inco- 
herent superposition of nonequilibrium and central ( "equilibrium" ) solutions 
of @ 

dN ch (y,t = T int ) »rl d / v 

j- = N C hRi(y, T int) 

+N 2 ch R 2 (y, r mt ) + N2R eq (y, r mt ). (8) 

The results for pseudorapidity distributions of produced charged hadrons 
in central Au + Au collisions at at two RHIC energies are shown in Fig. [3] in 
comparison with PHOBOS data [22l [23] . In the y 2 — minimization, the three- 
sources model yields excellent agreement with the data. Here the overall 
normalization is taken from the data, and the fit parameters are the time 
parameter (that determines the mean values < yi >2 >), the widths r 12 ,r eg , 
and the number of produced particles in the central source N^. 

At RHIC energies, the multiplicity density at midrapidity has still a sub- 
stantial contribution from the overlapping fragmentation sources. At 0.13 
TeV, the contribution from the three sources at rj = is about equal, at 0.2 
TeV the midrapidity source is larger (58%), but the fragmentation sources 
still contribute 21% each. 
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It should be mentioned that there exist detailed microscopic calculations 
of fragmentation sources from gq — >■ q and qg — >• q diagrams by Szczurek 
et al. [26j 127] for pion production in proton-proton, and nucleus-nucleus 
collisions at SPS and RHIC energies. These processes are also responsible 
for the observed differences [28] in the production of positively and negatively 
charged hadrons, in particular, pions. An extension of these calculations to 
LHC energies is very desirable. 

Within the 3-sources RDM, we had presented predictions at LHC energies 
of 5.52 TeV in [21] that were included in [31 Hj. The total number of produced 
charged hadrons had been extrapolated with log(y / sjv r /v) to obtain 26.5 * 
Np a rt at 5.52 TeV, with the number of participants N part . Based on this 
assumption, the calculated RDM-pseudorapidity distribution function turned 
out to underestimate the midrapidity result that is expected using the recent 
ALICE 2.76 TeV data point [2] by a factor of 2.7. 

I have now chosen to adjust the RDM parameters such that the ALICE 
midrapidity value at 2.76 TeV is reproduced, 1601±60 [2] (1584±4 (stat.)±76 
(sys.) in []]). On this basis, the RDM distribution functions at 2.76 and 5.52 
TeV can be calculated. 

With the extrapolation of the time parameter and the partial widths 
ri,2,eg from Fig. [2j plus corresponding extrapolations of the number of pro- 
duced particles in fragmentation and central sources as functions of log (-^/s aw) 
given in Table [TJ the results are shown in Fig. [3] The main uncertainty is in 
the extrapolation of the particle content of the fragmentation sources since 
the content of the central source is essentially fixed by the ALICE midrapidity 
data point. The calculation at 5.52 TeV is performed based on an extrapo- 
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lation of the multiplicity density at midrapidity with log (^s^n) that yields 
dN/drj ~ 1860 at midrapidity. 

At LHC energies, the overall scenario changes even more in favor of par- 
ticle production from the midrapidity source. The bulk of the midrapidity 
density is generated in the central source (93%), there is only a small overlap 
of the fragmentation sources at midrapidity as shown in Fig. |4} 

In a comparison with calculations at LHC energies that do not include 
the Jacobian transformation as displayed by the dotted curves in Fig. 3J4, it 
is evident that the midrapidity dip structure is essentially determined by the 
interplay of the three sources for particle production, and only marginally 
influenced by the transformation from y— to rj— space at these high energies. 
The central distribution including the Jacobian has no dip at LHC energies, 
but only a slight reduction in absolute magnitude at midrapidity, as shown 
by the dashed curve in Fig. |4j 

The smallness of the fragmentation sources at midrapidity is in qualitative 
agreement with results of a microscopic model that we had developed in [29] 
to investigate net-baryon distributions at LHC energies. In that approach, 
the net-baryon yield at large rapidities is calculated from the interaction of 
valence quarks with the gluon condensate in the respective other nucleus. Ex- 
tending the model to the midrapidity region |30j, a net-baryon midrapidity 
density dN/dy(y=0)~4 is obtained at 5.52 TeV, corresponding to a midra- 
pidity density of 12 valence quarks - as opposed to a total of 1248 valence 
quarks in the system. Hence the charged-hadron production from valence 
quark - gluon interactions at LHC energies can be expected to be very small 
in the midrapidity region. 
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4. Conclusion 

Based on the description of charged-hadron pseudorapidity distributions 
in central collisions of heavy symmetric systems at RHIC energies in a non- 
equilibrium-statistical model, I have presented predictions of pseudorapidity 
distributions of produced charged hadrons for central Pb + Pb collisions at 
LHC energies of 2.76 and 5.52 TeV. These rely on the extrapolation of the 
transport parameters in the relativistic diffusion model (RDM) with increas- 
ing center-of-mass energy. 

In a three-sources model, the midrapidity source that is associated with 
gluon-gluon collisions accounts for about 93% of the charged-particle mul- 
tiplicity density measured by ALICE at midrapidity in Pb + Pb collisions 
at 2.76 TeV. The fragmentation sources that correspond to particles that 
are mainly generated from valence quark - gluon interactions are centered 
at relatively large values of pseudorapidity (< 771,2 >~< 2/1,2 >— =F4.1) and 
hence, these contribute only marginally to the midrapidity yield. 

Since the Jacobian transformation from rapidity to pseudorapidity space 
is close to 1 at LHC energies due to the large mean transverse momenta, 
the size of the midrapidity-dip in the pseudorapidity distribution function is 
essentially determined by the relative particle content in the three sources, 
not by the Jacobian. Small corrections of the extrapolated values for the 
number of produced particles in the fragmentation sources may be required 
once the measured distributions become available from CMS, ATLAS and 
ALICE at both LHC energies. 
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